In silico reversal of repeat-induced point mutation (RIP) identifies the origins of repeat families and uncovers obscured duplicated genes

Background: Repeat-induced point mutation (RIP) is a fungal genome defence mechanism guarding against transposon invasion. RIP mutates the sequence of repeated DNA and over time renders the affected regions unrecognisable by similarity search tools such as BLAST. Results: DeRIP is a new software tool developed to predict the original sequence of a RIP-mutated region prior to the occurrence of RIP. In this study, we apply deRIP to the genome of the wheat pathogen Stagonospora nodorum SN15 and predict the origin of several previously uncharacterised classes of repetitive DNA. Conclusions: Five new classes of transposon repeats and four classes of endogenous gene repeats were identified after deRIP. The deRIP process is a new tool for fungal genomics that facilitates the identification and understanding of the role and origin of fungal repetitive DNA. DeRIP is open-source and is available as part of the RIPCAL suite at http://www.sourceforge.net/projects/ripcal

Gene validation and remodelling using proteogenomics of Phytophthora cinnamomi, the causal agent of Dieback

Phytophthora cinnamomi is a pathogenic oomycete that causes plant dieback disease across a range of natural ecosystems and in many agriculturally important crops on a global scale. An annotated draft genome sequence is publicly available (JGI Mycocosm) and suggests 26,131 gene models. In this study, soluble mycelial, extracellular (secretome), and zoospore proteins of P. cinnamomi were exploited to refine the genome by correcting gene annotations and discovering novel genes. By implementing the diverse set of sub-proteomes into a generated proteogenomics pipeline, we were able to improve the P. cinnamomi genome annotation. Liquid chromatography mass spectrometry was used to obtain high confidence peptides with spectral matching to both the annotated genome and a generated 6-frame translation. Two thousand seven hundred sixty-four annotations from the draft genome were confirmed by spectral matching. Using a proteogenomic pipeline, mass spectra were used to edit the P. cinnamomi genome and allowed identification of 23 new gene models and 60 edited gene features using high confidence peptides obtained by mass spectrometry, suggesting a rate of incorrect annotations of 3% of the detectable proteome. The novel features were further validated by total peptide support, alongside functional analysis including the use of Gene Ontology and functional domain identification. We demonstrated the use of spectral data in combination with our proteogenomics pipeline can be used to improve the genome annotation of important plant diseases and identify missed genes. This study presents the first use of spectral data to edit and manually annotate an oomycete pathogen

The Arabidopsis KH-Domain RNA-Binding Protein ESR1 Functions in Components of Jasmonate Signalling, Unlinking Growth Restraint and Resistance to Stress

Glutathione S-transferases (GSTs) play important roles in the protection of cells against toxins and oxidative damage where one Arabidopsis member, GSTF8, has become a commonly used marker gene for early stress and defense responses. A GSTF8 promoter fragment fused to the luciferase reporter gene was used in a forward genetic screen for Arabidopsis mutants with up-regulated GSTF8 promoter activity. This identified the esr1-1 (enhanced stress response 1) mutant which also conferred increased resistance to the fungal pathogen Fusarium oxysporum. Through positional cloning, the ESR1 gene was found to encode a KH-domain containing RNA-binding protein (At5g53060). Whole transcriptome sequencing of esr1-1 identified altered expression of genes involved in responses to biotic and abiotic stimuli, hormone signaling pathways and developmental processes. In particular was an overall significant enrichment for jasmonic acid (JA) mediated processes in the esr1-1 down-regulated dataset. A subset of these genes were tested for MeJA inducibility and we found the expression of some but not all were reduced in esr1-1. The esr1-1 mutant was not impaired in other aspects of JA-signalling such as JA- sensitivity or development, suggesting ESR1 functions in specific components of the JA-signaling pathway. Examination of salicylic acid (SA) regulated marker genes in esr1-1 showed no increase in basal or SA induced expression suggesting repression of JA-regulated genes is not due to antagonistic SA-JA crosstalk. These results define new roles for KH-domain containing proteins with ESR1 unlinking JA-mediated growth and defense responses.http://doi.org/10.1371/journal.pone.0126978</a

“CATAStrophy,” a Genome-Informed Trophic Classification of Filamentous Plant Pathogens – How Many Different Types of Filamentous Plant Pathogens Are There?

The traditional classification of fungal and oomycete phytopathogens into three classes – biotrophs, hemibiotrophs, or necrotrophs – is unsustainable. This study highlights multiple phytopathogen species for which these labels have been inappropriately applied. We propose a novel and reproducible classification based solely on genome-derived analysis of carbohydrate-active enzyme (CAZyme) gene content called CAZyme-Assisted Training And Sorting of -trophy (CATAStrophy). CATAStrophy defines four major divisions for species associated with living plants. These are monomertrophs (Mo) (corresponding to biotrophs), polymertrophs (P) (corresponding to necrotrophs), mesotrophs (Me) (corresponding to hemibiotrophs), and vasculartrophs (including species commonly described as wilts, rots, or anthracnoses). The Mo class encompasses symbiont, haustorial, and non-haustorial species. Me are divided into the subclasses intracellular and extracellular Me, and the P into broad and narrow host sub-classes. This gives a total of seven discrete plant-pathogenic classes. The classification provides insight into the properties of these species and offers a facile route to develop control measures for newly recognized diseases. Software for CATAStrophy is available online at https://github.com/ccdmb/catastrophy. We present the CATAStrophy method for the prediction of trophic phenotypes based on CAZyme gene content, as a complementary method to the traditional tripartite “biotroph–hemibiotroph–necrotroph” classifications that may encourage renewed investigation and revision within the fungal biology community.</p

Dwarf Cepheids in the Carina Dwarf Spheroidal Galaxy

We have discovered 20 dwarf Cepheids (DC) in the Carina dSph galaxy from the analysis of individual CCD images obtained for a deep photometric study of the system. These short-period pulsating variable stars are by far the most distant (~100 kpc) and faintest (V ~ 23.0) DCs known. The Carina DCs obey a well-defined period-luminosity relation, allowing us to readily distinguish between overtone and fundamental pulsators in nearly every case. Unlike RR Lyr stars, the pulsation mode turns out to be uncorrelated with light-curve shape, nor do the overtone pulsators tend towards shorter periods compared to the fundamental pulsators. Using the period-luminosity (PL) relations from Nemec et al. (1994 AJ, 108, 222) and McNamara (1995, AJ, 109, 1751), we derive (m-M)_0 = 20.06 +/- 0.12, for E(B-V) = 0.025 and [Fe/H] = -2.0, in good agreement with recent, independent estimates of the distance/reddening of Carina. The error reflects the uncertainties in the DC distance scale, and in the metallicity and reddening of Carina. The frequency of DCs among upper main sequence stars in Carina is approximately 3%. The ratio of dwarf Cepheids to RR Lyr stars in Carina is 0.13 +/- 0.10, though this result is highly sensitive to the star-formation history of Carina and the evolution of the Horizontal Branch. We discuss how DCs may be useful to search effectively for substructure in the Galactic halo out to Galactocentric distances of ~100 kpc.Comment: 20 pages of text, 7 figure

OcculterCut: A comprehensive survey of AT-rich regions in fungal genomes.

We present a novel method to measure the local GC-content bias in genomes and a survey of published fungal species. The method, enacted as "OcculterCut" (https://sourceforge.net/projects/occultercut), identified species containing distinct AT-rich regions. In most fungal taxa, AT-rich regions are a signature of repeat-induced point mutation (RIP), which targets repetitive DNA and decreases GC-content though the conversion of cytosine to thymine bases. RIP has in turn been identified as a driver of fungal genome evolution, as RIP mutations can also occur in single-copy genes neighbouring repeat-rich regions. Over time RIP perpetuates 'two speeds' of gene evolution in the GC-equilibrated and AT-rich regions of fungal genomes. In this study, genomes showing evidence of this process are found to be common, particularly among the Pezizomycotina. Further analysis highlighted differences in amino acid composition and putative functions of genes from these regions, supporting the hypothesis that these regions play an important role in fungal evolution. OcculterCut can also be used to identify genes undergoing RIP-assisted diversifying selection, such as small, secreted effector proteins that mediate host-microbe disease interactions

Proteomic analysis of Rhizoctonia solani identifies infection-specific, redox associated proteins and insight into adaptation to different plant hosts

Rhizoctonia solani is an important root infecting pathogen of a range of food staples worldwide including wheat, rice, maize, soybean, potato and others. Conventional resistance breeding strategies are hindered by the absence of tractable genetic resistance in any crop host. Understanding the biology and pathogenicity mechanisms of this fungus is important for addressing these disease issues, however, little is known about how R. solani causes disease. This study capitalises on recent genomic studies by applying mass spectrometry based proteomics to identify soluble, membrane-bound and culture filtrate proteins produced under wheat infection and vegetative growth conditions. Many of the proteins found in the culture filtrate had predicted functions relating to modification of the plant cell wall, a major activity required for pathogenesis on the plant host, including a number found only under infection conditions. Other infection related proteins included a high proportion of proteins with redox associated functions and many novel proteins without functional classification. The majority of infection only proteins tested were confirmed to show transcript up-regulation during infection including a thaumatin which increased susceptibility to R. solani when expressed in Nicotiana benthamiana. In addition, analysis of expression during infection of different plant hosts highlighted how the infection strategy of this broad host range pathogen can be adapted to the particular host being encountered. Data are available via ProteomeXchange with identifier PXD002806

Deep proteogenomics; high throughput gene validation by multidimensional liquid chromatography and mass spectrometry of proteins from the fungal wheat pathogen Stagonospora nodorum

Background Stagonospora nodorum, a fungal ascomycete in the class dothideomycetes, is a damaging pathogen of wheat. It is a model for necrotrophic fungi that cause necrotic symptoms via the interaction of multiple effector proteins with cultivar-specific receptors. A draft genome sequence and annotation was published in 2007. A second-pass gene prediction using a training set of 795 fully EST-supported genes predicted a total of 10762 version 2 nuclear-encoded genes, with an additional 5354 less reliable version 1 genes also retained. Results In this study, we subjected soluble mycelial proteins to proteolysis followed by 2D LC MALDI-MS/MS. Comparison of the detected peptides with the gene models validated 2134 genes. 62% of these genes (1324) were not supported by prior EST evidence. Of the 2134 validated genes, all but 188 were version 2 annotations. Statistical analysis of the validated gene models revealed a preponderance of cytoplasmic and nuclear localised proteins, and proteins with intracellular-associated GO terms. These statistical associations are consistent with the source of the peptides used in the study. Comparison with a 6-frame translation of the S. nodorum genome assembly confirmed 905 existing gene annotations (including 119 not previously confirmed) and provided evidence supporting 144 genes with coding exon frameshift modifications, 604 genes with extensions of coding exons into annotated introns or untranslated regions (UTRs), 3 new gene annotations which were supported by tblastn to NR, and 44 potential new genes residing within un-assembled regions of the genome. Conclusion We conclude that 2D LC MALDI-MS/MS is a powerful, rapid and economical tool to aid in the annotation of fungal genomic assemblies

Mass-spectrometry data for Rhizoctonia solani proteins produced during infection of wheat and vegetative growth

© 2016. Rhizoctonia solani is an important root infecting pathogen of a range of food staples worldwide including wheat, rice, maize, soybean, potato, legumes and others. Conventional resistance breeding strategies are hindered by the absence of tractable genetic resistance in any crop host. Understanding the biology and pathogenicity mechanisms of this fungus is important for addressing these disease issues, however, little is known about how R. solani causes disease. The data described in this article is derived from applying mass spectrometry based proteomics to identify soluble, membrane-bound and culture filtrate proteins produced under wheat infection and vegetative growth conditions. Comparisons of the data for sample types in this set will be useful to identify metabolic pathway changes as the fungus switches from saprophytic to a pathogenic lifestyle or pathogenicity related proteins contributing to the ability to cause disease on wheat. The data set is deposited in the PRIDE archive under identifier PRIDE: PXD002806

A first genome assembly of the barley fungal pathogen Pyrenophora teres f. teres

Background: Pyrenophora teres f. teres is a necrotrophic fungal pathogen and the cause of one of barley’s most important diseases, net form of net blotch. Here we report the first genome assembly for this species based solely on short Solexa sequencing reads of isolate 0-1. The assembly was validated by comparison to BAC sequences, ESTs, orthologous genes and by PCR, and complemented by cytogenetic karyotyping and the first genome-wide genetic map for P. teres f. teres. Results: The total assembly was 41.95 Mbp and contains 11,799 gene models of 50 amino acids or more. Comparison against two sequenced BACs showed that complex regions with a high GC content assembled effectively. Electrophoretic karyotyping showed distinct chromosomal polymorphisms between isolates 0-1 and 15A, and cytological karyotyping confirmed the presence of at least nine chromosomes. The genetic map spans 2477.7 cM and is composed of 243 markers in 25 linkage groups, and incorporates SSR markers developed from the assembly. Among predicted genes, non-ribosomal peptide synthetases and efflux pumps in particular appear to have undergone a P. teres f. teres-specific expansion of non-orthologous gene families. Conclusions: This study demonstrates that paired-end Solexa sequencing can successfully capture coding regions of a filamentous fungal genome. The assembly contains a plethora of predicted genes that have been implicated in a necrotrophic lifestyle and pathogenicity and presents a significant resource for examining the bases for P. teres f. teres pathogenicity